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ABSTRACT 


A  theory  of  Fabry-Perot  resonances  designed  to  be  useful 
at  small  Fresnel  numbers  is  given.  This  *w,eory  is  applied  to  the 
case  where  the  interferometer  mirrors  are  or  arbitrary  curvature 
and  have  circular  apertures.  Asymptotic  formulas  are  derived 
for  the  diffraction  output  and  phase  shifts  in  the  limit  of  small 
Fresnel  numbers.  These  formulas  demonstrate  the  mode  discrimina¬ 
tion  properties  of  the  interferometers  when  operated  in  this  re¬ 
gime.  Numerical  results  are  presented  covering  a  wider  range  of 
Fresnel  numbers. 


11 . 


1. 


INTRODUCTION 

The  resonance  properties  of  Fabry-Perot  (F.P.)  inter¬ 
ferometers  at  Fresnel  numbers  of  the  order  of,  and  greater  than, 
unity  have  been  treated  in  a  number  of  papers, The  recent 
use  of  "diffraction  output  coupling"  as  an  efficient  method  of 
transverse  mode  suppression  when  the  Fresnel  numbers  are  less 
than  unity^  has  motivated  us  to  extend  the  theory  so  as  to 
describe  resonators  operating  in  this  range.  We  have  previously 
presented  numerical  results  for  small  Fresnel  numbers,  on  diffrac¬ 
tion  outputs  and  phase  shifts  for  F.P.  interferometers  consisting 
of  plane  mirrors  with  circular  apertures. ^ 

In  this  article  we  give  a  theory  of  F.P.  resonances 
designed  to  be  useful  at  small  Fresnel  numbers.  This  theory 
is  applied  to  the  case  where  the  interferometer  mirrors  are  of 
arbitrary  (small)  curvature  and  have  circular  apertures. 

Asymptotic  formulas  are  derived  for  the  diffraction  output  and 
phase  shifts  in  the  limit  of  small  Fresnel  numbers.  These  formulas 
demonstrate  the  mode  discrimination  properties  of  the  interferom¬ 
eters  when  operated  in  this  regime.  Numerical  results,  in  addition 
to  those  previously  reported v  '  ,  are  presented  on  the  diffraction 
output  and  phase  shifts  for  curved  circular  mirrors  for  small 
Fresnel  numbers.  The  theory  is  presented  in  such  a  way  (in  terms 
of  a  canonical  integral  equation)  that  it  applies  to  a  general 
class  of  aperture  shapes. 


FORMULATION  OF  THE  PROBLEM 


The  geometry  which  is  considered  consists  of  a  symmetric 
pair  of  mirrors,  each  of  radius  of  curvature,  R,  whose  subtending 
planes  are  separated  along  the  axis  of  the  resonator  by  a  distance 
L,  The  aperture  shapes  are  taken  to  be  circles  of  radius  a  .  (See 
Figure  1).  We  are  interested  in  the  resonance  properties  of  this 
electromagnetic  cavity,  in  particular,  the  complex  resonant  freq- 
uencies  and  the  corresponding  field  configurations.  The  informa¬ 
tion  contained  in  the  complex  resonant  frequencies  can  be  expressed 
in  terms  of  the  single  pass  diffraction  output  and  phase  shift  of 
the  resonant  fields. 

Consider  an  arbitrary  field  configuration  at  plane  (1) 
of  Figure  1.  This  field  may  be  regarded  as  the  sum  of  two  fields, 
one  moving  to  the  right,  ^f+(x,y),  and  one  to  the  left  ^f_(x,y). 
The  field  moving  to  the  right  reflects  from  the  other  mirror  and 
results  in  a  configuration  v  /f_(x,y),  representing  a  field  moving 
to  the  left  at  plane  (2) .  The  resonances  of  this  open  structure 
are  defined  by  the  requirement  that  after  a  double  transit  the 
field  reproduces  itself,  or,  after  a  single  pass  it  reproduces 
itself  up  to  the  sign  change  -1,  that  is 


(2^f_(x,y)  =  -  ^f _(x,y) 

More  generally,  we  have: 

^f_(x,y)  =  relK(k)L  (1)f^(x,y), 


(1) 

(2) 


where  K(k)L  is  a  complex  number,  independent  of  x  and  y,  but  a 
function  of  the  geometry  and  propagation  constant  and  r  is  the 
mirror  reflectivity.  The  single  pass  phase  shift,  tf,  is  given 
by : 

- 

The  method  of  analysis  can  be  extended  to  other  aperture  shapes, 
m  particular,  rectangles. 

Because  the  structure  is  open,  the  _esonant  frequencies  are 
complex.  The  imaginary  part  of  the  frequency  is  the  inverse  of 
the  lifetime  of  the  corresponding  f,mode". 


3. 


0  -  Re  K(k) L ,  (3) 

and  the  single  pass  fractional  power  output  (due  to  diffraction) , 

2 

1-  j~yj  ,  is  determined  by 


(4) 


The  resonant  values  of  the  frequency.  Re  ck,  and  gain 
coefficient,  Imk,  are  determined  by  the  resonance  condition  (1) , 
or 


re«(k)L  .  tl 

The  quantity  of  interest,  K(k) ,  is  determined  by  the 
solution  of  the  integral  equation  which  is  obtained  by  propagating 
the  field  from  plane  (1)  to  plane  (2)  and  using  equation  (2) . 

The  following  approximations  allow  a  relatively  simple  formulation 
of  the  problem: 

(1)  The  mirror  radius  of  curvature,  R,  is  large  enough  that 
it  is  consistent  to  take  the  electric  (or  magnetic) 
field  vector  to  be  always  perpendicular  to  the  resonator 
axis.  This  reduces  the  real  situation  to  a  scalar 
problem  where  the  total  field  is  determined  by  a  single 
transverse  component. 

(2)  For  a  field  on  a  plane  intersecting  the  rim  of  a  mirror 
(such  as  plane  (1)  or  (2)  in  Figure  1)  moving  toward  the 
mirror,  the  reflected  field  moving  away  from  the  mirror, 
at  the  point  (x,y)  on  this  plane  is: 

2  2  2 

(a)  zero,  for  x  -  y  >  a  (Kirchoff-approximation, 

(b)  equal  in  amplitude  to  the  reflection  coefficient, 
r,  times  the  incident  field,  and  increased  in 

f  a2  »  x2^  2j~ 

phase  factor  by  the  amount  exp  !2  i  jl-  — 7^—, 

L  '  'j 


A. 


O  *>  O 

for  <  a\  This  phase  is  the  amount  that 

would  be  accumulated  by  a  ray  moving  parallel  to 
the  axis. 

9  2 

(c.)  (L/a)  »  (a  /L%)  =  Fresnel  number. 

Because  of  the  axial  symmetry  the  solutions  must  be  of  the  form: 

f(x,y)  =  Rn(p)e^n®,  where 
x  =  ap  cos  0 
y  =  ap  sin  9. 

Using  assumptions  (a) ,  (b)  and  (c)  and  the  resonance  condition 

■ k 

(1) ,  we  obtain  the  integral  equation: 


7 


r  m(p>  = p'dp' 


n,m  ”  n,m 


R  (p  ’ )  e 
n,mVh  ' 


i-Np2  ei-N ; 1 -2L/r: p ' 2Jn(2-Nf  p ' )  (5) 


The  propagation  factors  exp(iKn  ^L)  are  obtained  from  the  eigen¬ 
values  't  of  equation  5  by  the  relation: 

exp(iK  L)  =  ,  2~N  exp’ikL  -  i7(n-l) . .  (6) 

r  n,m  nsm  / 

In  the  above  relations,  the  Fresnel  number,  N,  is  given  by: 

N  -  a2/!/'  (7) 


Canonical  Integration  Equation 

The  integral  equation  (5)  may  be  identified  with  a 
canonical  form: 

b 

r-_  (x)  =  '  0(x)  K(-,xx’)  0(x’)  .(x’)  (x1)  dx’ .  (8) 

a 

This  same  form  arises  when  the  aperture  shape  corresponds  to 
separable  coordinate  systems  other  than  circular,  in  particular, 
rectangular.  For  the  specific  case  of  circular  apertures  the 

•k 

Set  Appendix  t\. 


5. 


particular  identifications  that  are  necessary  in  order  to  convert 
the  integral  equation  (5)  into  the  form  (8)  are: 

X  =  p2,  (9a) 


K(  =  )  =  n!  r"/2  Jn(2  -;) 

,  s2  (9b) 

-  1  '  TT'fn-'n  +  2'.(n*lhn+2)  . 

n 

“T 

t (x)  =  x  Rn(x)  exp: -iTNLx/Rj ,  (9c) 


•'  -  2n!  Vn/('N>n’ 

(9d) 

0(x)  =  exp. i~Nx(l-L/R) ] , 

(9e) 

(x)  =  xn, 

(9f) 

-  (~N)2, 

(9g) 

a  =  0,  b  =  1 . 

(9h) 

Two  methods  for  obtaining  solutions  to  equation  (8)  will  be 
presented.  One  is  particularly  convenient  for  getting  the  asymp¬ 
totic  behavior  in  the  limit  of  small  ~  and  expansions  in  power 
series  of  the  Fresnel  number  useful  for  ~N  <  1 .  The  other  is 
useful  as  a  basis  for  machine  computations  in  order  to  obtain 
the  resonance  properties  for  values  of  the  Fresnel  number  up 
to  unity. 


4. 


6. 


ASYMPTOTIC  SOLUTIONS 


In  this  section  we  investigate  the  eigenvalues,  F,  and 
the  eigenfunctions,  (x) ,  of  the  canonical  integral  equation  (8) 
in  the  limit  q  0.  Correspondence  will  be  made  with  the  asymp¬ 
totic  properties  of  the  resonances  of  the  circular  aperture  F.P. 
interferometer  in  the  limit  N  -*  0,  by  means  of  the  relations  (9a) 
to  (9h)  . 

The  eigenvalues  F,  eigenfunctions  i  (x)  and  the  kemal 
K(qxx')  are  expanded  in  power  series  in  the  parameter  q: 


m 


m=o 


m 


(10) 


.  (x>  =  E  (11) 

m-o 

and 

K(-,xx')  ==  r.  k  xm  x’m  1m.  (12) 

m=o 


When  these  three  expansions  are  substituted  into  the  integral 
equation  (8),  we  obtain,  upon  equating  coefficients  of  equal 
powers  of  q: 


s 

z  'K 


'm(x)  7s-m  -  0(x)  "  kx,  ^  Cs.m* 

m=o  m=o 


where 


and 


Cs  m  =  dx 1  .(x')  x'm  0(x’ )  .(::’), 

S  -  TP. 


s  -=  0,  1 ,  2 ,  ... 


(13a) 


(13b) 


(13c) 


We  observe,  from  equations  (13)  that  either  (x)  is  of  the  form: 


7. 


m 


(x)  =  0(x) 


(Polynomial  of  degree  m) , 


or  F  =0.  Calling  the  first  solution  ,  and  the  second 

m  °  n\  m(x) 

'  '•  ,  we  find,  either  ' 


m(x) 


m(x) 


is  of  the  form: 


(1) 


m 


(x)  ~  *  (Polynomial  of  degree  m  1) , 


or  T,  =0.  Continuing  in  this  way  we  see  that  there  are  an  in- 
1  /  *  \ 


finite  number  of  modes 

(i) 


a). 


(x) 


,  with  the  corresponding  eigen¬ 


values  v  'F,  which  have  the  expansions: 


(«, 


'(x) 


X 

-  V 

—  t— 

m=o 


(i). 


’m(x) 

(i) 


m=o 


(i)  ,  m 

Vm(x)  ^  ’ 

(14) 

^  (i)P  x' 

^  p_  ^  X  i 

f)  in  •  4 

(15) 

i-O  9 

(i)-  m+i 

'  m  * 

(16) 

Substituting  these  relations  into  the  integral  equation  (8)  and 
equating  coefficients  of  equal  powers  of  tj,  and  then  equating 
coefficients  of  equal  powers  of  x,  we  finally  obtain  the  relations: 


where 


s-m 


n=o 


i-n,m 


s-m-*-i 

k  I  u)d  t 

m  rs-m,n  n-f-m’ 

n-o  * 


(17) 


L 

t  = 
m  a 

02(x) 

xm  —  (x)  dx, 

(18) 

^n,m 

0  for 

n  negative 

(19) 

s  - 

0,  1, 

2,  ... 

(20) 

m  -  0 , 

1,  2, 

. . .  s-1,  s. 

(21) 

i* 


*>  4  /-n«^  o  /17  \  /Ol  \  ^ _ i  ^ 

4-.  vjmu  V  x  v^*io  \x/  /  Uli.  wu^u^a;  v  Uinp  x  v  u 


£:  1  ••  ^  ^  *-l>  r»  r  /->  1  1 1  f  i'  rtn  c  c  1  n  p  o 

c  x  y  ut  u c  t  utn  iv-  tu  v-  «ju  i.  u  l.  t  vu^  ^  ttv  w 


they  determine  the  quantities  *  'p  and  '  '  which  through 

m ,  /  m  ° 

equations  (14),  (15)  and  (16)  give  the  eigenfunctions  and  eigen¬ 
values.  The  quantities  t  and  k  are  known.  The  wave  functions 

m  m 

are  determined  up  to  an  arbitrary  normalization  factor.  A  con¬ 
venient  normalization,  which  Is  used  below,  is: 


-O)  -  0(0). 


(22) 


This  normalization  implies: 


(23) 


The  coefficients  necessary  in  order  to  get  the  first 
two  terms  in  the  expansions  of  ^  -  >.  and  are,  for  the 

\X) 

cases  i  -  0,  and  i  -  1: 

CASE  1  j  ^  0 

°r  --  k  t 
o  o  o 

°r  -  k-,t,2/t_ 

i  i  i  o 


r>  ,  r\ 


?Ur?oo  Rlcl/koto 


CASE  2  i  =  1 


1_  ki  ,  9 

-  — -  (t0t  -  t, ") 
o  t  2  o  i 


(24a) 

(24b) 


(24c) 


(25a) 


(lp01/lpoo>  -  -  to/t! 


(25b) 


^lpi2/lpoo)  =  Pgtj  (t2Cl  '  t3to) 


(25c) 


(lpn/lpoo>  s  trr- 

i  o 

"1P  1  p 

u  __u  t  ,  1112 

kl  ip  r2  lPni 

L  r'l  01 


P12!  C2 


\  Poo/  ^ 


(25d) 


(2  5t) 


9. 


We  now  specialize  to  the  case  of  circular  apertures. 

Using  the  relations  (9a)  -  (9h),  and  equation  (6), we  have,  for  the 

propagation  factor: 


exp  i  K  L 
l  n,m 


/  x[vn+2m4-l 

n:  '  o 


exp 


ikL-iJ(n^l) 


2-iNL/R 


(26) 


The  quantity  mro  is  determined  by  equation  (17)  as  a  function  of 
the  lc‘s  and  t‘s.  The  latter  quantities  are  given,  asymptotically, 
by: 


1 

n-j+1 


exp 


i2-(l-L/R; 


0(N2) 


(27) 


The  former  are  given  by: 


k j  =  (- 3 )  ^  n ! /  (j  '.)  (n  -  j) 


(28) 


As  the  Fresnel  number  gets  very  small,  the  quantity  approaches 
a  constant  (independent  of  N)  and  hence  the  dependence  on  the  Fresnel 
number  of  the  single  pass  diffraction  loss  is  determined  by: 


|  exp 
I 


r 


K  L 
n,m 


N-*0 


N 


2n*4m-2 


(29) 


This  relation  illustrates  the  strong  mode  discrimination  in  the 
limit  of  small  Fresnel  numbers. 


4. 


10. 


MORE  GENERAL  EXPANSIONS 


The  complete  expansion  in  powers  of  the  Fresnel  number 
for  the  lowest  modes  (m=0)  for  any  value  of  n  can  be  obtained 
from  sample  recurrence  relations.  ITius,  we  have  from  equation 
(17)  for  i  =  0,  m  -  0 


s 


s 


n=o 


s ,  n 


(30) 


where 


(31) 


and  the  quantities  p  for  n  ^  s  are  related  to  previously 

s  2  n 

determined  quantities  by  the  expression: 


P 


s,n 


s-n 


J^o 


s-n,  j 


1  s-n 

1»1  J  Ps-J,n, 


(32) 


which  is  obtained  from  equation  17  for  i=0,  by  setting  m^n. 

The  con*e spending  complete  expansion  for  the  higher 
modes  (i  ">  0)  cannot  be  obtained  in  terms  of  simple  recurrence 
relations,  but  rather  require  the  solution  of  simultaneous  algebraic 
equations.  In  order  to  obtain  numerical  results  for  the  case 
that  the  Fresnel  number  is  not  much  less  than  unity  an  alternative 
method  of  analysis  to  that  described  above  may  be  employed.  This 
method  treats  groups  of  modes  on  an  equal  basis  and  hence  while 
not  as  satisfactory  for  the  lowest  mode  is  more  easily  mechanized 
for  machine  computation. 

Referring  to  the  canonical  integral  equation  (8),  for 
the  Mth  degree  of  approximation  we  represent  the  kernal  K(  ),  by 
the  Mth  degree  polynomial: 

M 


n  o 


(33) 


11 . 


The  integral  equation  becomes: 

7  .  (x)  -  :  knxrSn  0(x)  0(x‘)  .(x’)  (x')dx'  . 

n-o  a 

The  solutions  of  this  integral  equation  must  be  of  the  form: 

M 

(x)  =  0(x)  I  P  Xn  . 
n=o 


(34) 


(35) 


The  coefficients  are  determined  by  substituting  (35)  into  (34) 
and  equating  coefficients  of  equal  powers  of  x.  In  this  way  one 
gets  the  M^l  homogeneous  equations  for  the  M--1  unknowns  P  ,  •••  P^: 


where,  as  before: 


M 


m=o 


"n-m  m’ 


t 


n 


b  9 

'  0z(x)  xn  . (x) dx . 

a 


(36) 


(37) 


The  M-l  eigenvalues  7,  are  solutions  of  the  determintal  equation: 


t 

n-m 


0 


(38) 


This  method  of  solution  has 


been  described  by  She  and  Heffner 


12. 


RESULTS 


From  equations  (24)  to  (28)  for  the  asymptotic  propaga¬ 
tion  factors  for  the  first  four  modes  are  given  by: 

exP  - 1  K0#o  L]  =  (~N)  exp  i;<J  -  t  "  ~N  g- 

exp.i  Q  L"  =  ^-(^N)^  exp  i'^  -  -  -  ^  ~Ng' 

exp.i  K0jl  L'.  -  -^-("N)3  ex?  1  5  ’  J  *  "Mg' 

exp!  i  K|  x  L]  =  -^(“N)4  exp  i^  -  y|  ~Ng  ' * 

where  ^  is  the  geometrical  phase  shift  given  by 

$  -  kL  *  2-NL/R, 

and 

g  -  1  -  L/R. 

It  is  generally  txrue  that  the  phase  shift,  relative  to  the  geo¬ 
metrical  phase  shift  is  linear  in  the  quantity  Ng  asymptotically 
for  small  enough  N.  As  was  mention  above,  the  asymptotic  depend¬ 
ence  on  the  Fresnel  number  of  the  single  pass  diffraction  loss  is 
determined  by: 


(39a) 


(39b) 


(39c) 


(39d) 


(39e) 


(39f) 


exp  i  K  L 
r  n,m 


t2n-4m  -2 


N-*o 


(40) 


The  results  for  a  wider  range  of  Fresnel  numbers  is 
illustrated  in  figures  (2)  to  (5) .  In  figures  (2)  and  (3)  we 

2  *  2 

have  plotted  hj  =  exp.i  KLl!^  for  the  first  two  modes,  (n^O,  m^O) 


(n=l ,  m=0)  respectively,  versus  the  Fresnel  number  a  /L  -.  The 

2 

reciprocal  of  J~,  j  is  the  single  pass  gain  required  to  reach 
the  threshold  of  oscillation  in  the  mode.  The  fractional  dif- 
fracted  power  output  per  transit  is  1  -  j,j^.  The  quantity 
g  -  (1  -  L/R)  measures  the  mirror  curvature.  The  case  g  --  0 


13. 


corresponds  to  a  confocal  arrangement  and  the  case  g  =  1  to  flat 
mirrors.  For  small  enough  values  of  the  Fresnel  number  the  losses 
are  insensitive  to  the  value  of  g.  The  detailed  dependence  on  g 
for  the  larger  values  of  N  is  given  in  expanded  scale  in  the  lower 

right  hand  corners  of  figures  (2)  and  (3) .  It  is  seen,  that  for 

\ 

the  losses,  the  asymptotic  expressions  (39)  become  quite  accurate 
for  Fresnel  numbers  less  than  0.1  for  the  lowest  mode  and  0.2  for 
the  first  excited  mode. 

In  figures  (4)  and  (5)  we  have  plotted  the  phase  shift 
(relative  to  the  geometrical  phase  shift  versus  the  Fresnel  number 
for  the  first  two  modes.  It  is  seen  that  the  asymptotic  expres¬ 
sions  (39)  give  the  phase  shift  quite  accurately  for  Fresnel  numbers 
less  than  0.1. 
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APPENDIX  A 


DERIVATION  OF  RESONANCE  INTEGRAL  EQUATION 

If  a  field  distribution  is  given  over  the  plane  z'=0, 
then  the  field  at  the  point  (x,ytL) ,  which  satisfies  the  scalar 
Helmoltz  equation  is: 


f(x,y,L)  =  ,  dx'dy'  f(x',y',z')  G(x,y ,L ,x' ,y '  ,z’ ) 


(A-l) 


where  G(x,y,L;  x',y',z')  is  the  Greens  function  for  the  scalar 
Helmoltz  equation  which  vanishes  in  the  plane  z'=0  and  represents 
outgoing  waves.  This  Greens  function  is  given  by: 


1  [  eik  J(x-x* ) 2+(y-y ' ) 2+(L-z ' ) 2 
G(x,y,L,  x',y',z')  =  -  t~ 

\/(x~x')2+<y-y')2-(L-z')2 


ik  v/(x-x')2+(y-y')2-'-(L+z’)2 
V(x-x' )  2+(y-y' )  2+(L-z')  2 


(A-2) 


Thus,  equation  (1)  may  be  written 


-1  e1** 

f(x,y,L)  =  ikL  .dx'dy'  f(x',y',0)  — 7-  , 

L  R 


(A- 3) 


where  R 


(x-x')2  +  (y-y')2  +  L2. 


If  the  mirror  diameter,  2a  is  small  compared  to  the 
distance  L,  then  we  may  make  the  expansion: 


16. 


1  +  1  Lx-x,)"Hy«yt) 

7  7 

L" 


»  \  *• 


<  /  i 

o|<| 

\L  J 


In  the  argument  of  the  exponential  term,  the  dropping 


of  the  remainder  0  -^rj  is  justified  if 


kL  0 


L5 


«  " 


This  implies  the  requirement: 


N 


(A-4) 


To  get  the  form  of  equation  (5)  we  must  include  the 
additional  phase  shift  and  the  reflectivity  of  the  mirrors. 

Using  assumption  (2)  of  the  section  "Formulation  of  the  Problem" 
we  have 


(7\  _,ikL-i'  /2  2  1  2-  /n 

f  (r  -  if® _ *L_  WQ'  (1) 


f+(p,e)  «  -r  77 ^ 

[L+az/R]  > 


-■  f'dp;  d0’  ll'£.(;',#')exp(-2  i|f'2) 


exp  i 


[L-a^/R] 


f2  -  f ,2-2, f '  cos  (e-e')>. 


(A-5) 


Using  the  relation: 


2-(-i)n  Jn(z) 


,  inu  iz  , . 

-  .  du  e  e  cos  a,  (A-6) 


and  taking  N  «  R/ ' ,  we  get  the  form  of  equation  (5). 
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